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Motivated by a nice result shown by E. Linder, a detailed discussion on the choice of the expansion 
parameters of the Maclaurin series for the equation of states of a perfect fluid is presented in this 
paper. We show that their nice recent result is in fact a linear approximation to the full Maclaurin 
series as a power series of the parameter y = z/(l + z). The power series for the energy density 
function, the Hubble parameter and related physical quantities of interest are also shown in this 
paper. The method presented here will have significant application in the precision distance-redshift 
observations aiming to map out the recent expansion history of the universe. In addition, a complete 
analysis of all known advantageous parameterizations for the equation of states to high redshift is 
also presented. 



I. INTRODUCTION 

Recently there are great advances in our abilities in cosmological observations for the quest of exploring the expansion 
history of the universe. It carries cosmology well beyond determining the present dimensionless density of matter 
VL m and the present deceleration parameter go of Sandage [1]. One seeks to reconstruct the entire function a(t) 
representing the expansion history of the universe. Earlier on, cosmologists sought only a local measurement of the 
first two derivatives of the scale factor a, evaluated at a single time t$. One tries instead to map out the function 
determining the global dynamics of the universe in the near future. While many qualitative elements of cosmology 
follow merely from the form of the metric, (see Weinberg [2]), deeper understanding of our universe requires knowledge 
of the qualitative dynamics of the scale factor a(t). This echoes the transition of energy between components signifying 
the epoch of radiation domination to that of matter domination. This is also a key element in the growth of density 
perturbations into large structure. Yet until recently the literature only considered the Hubble constant and the 
deceleration parameter measured today. There are now a great numbers of cosmological observational tests that will 
be able to probe the function a(t) more completely throughout all ages of the universe (see Sandage [3] , Linder [4] , 
[5], Tegmark [6]). 

What one required is a probe capable of both precise and accurate enough observations. Indeed, a number of 
promising methods are being developed including the magnitude-redshift relation of Type la supernovae. The goal 
of mapping out the recent expansion history of the universe is well motivated. The thermal history of the universe, 
extending back through structure formation, matter-radiation decoupling, radiation thermalization, primordial nu- 
cleosynthesis, etc. is very important in the study of cosmology and particle physics, high energy physics, neutrino 
physics, gravitational physics, nuclear physics, and so on (see, e.g., Kolb and Turner [7]). 

Recent expansion history of the universe is similarly a very promising research focus with the discovery of the 
current acceleration of the expansion of the universe. This includes the study of the role of high energy field theories 
in the form of possible quintessence, scalar-tensor gravitation, higher dimension theories, brane worlds, etc in the very 
recent universe. The accelerated expansion is also important to the possible fate of the universe [8], [9], [10]. It is 
then important to consider the use of supernova observations to obtain the magnitude-redshift law back to z ~ 1.7 
and how to relate this to the scale factor-time behavior a(t) with the proposed Supernova/Acceleration Probe mission 
[11]. 

Therefore, the study of modelling different equation of state (EOS) derived from different theories plays an important 
role in the study of the very recent expansion history of our universe. Interestingly, there is a nice paper showing a 
new model that is capable of extracting important physics with the help of the linear parameterization in powers of an 
appropriate parameter z/(l + z) expanding the function uj of the equation of state [12], [13], [14], [15]. We will show 
in section II that one should treat the proposed new model as the leading term in the Maclaurin series associated 
with the EOS (cf. [14] and Fig. 1 of [15]). In section III, we will present the first four leading terms in computing the 
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energy density function p, the Hubble parameter H, and the conformal time 77 with the help of the Maclaurin series 
of the EOS. 



II. TAYLOR EXPANSION OF THE EQUATION OF STATE 

The equation of state is given by the relation ui = p/p. In addition, the field equations of the universe can be shown 
to be: 

d]np(z) = 3(u> + l)d]n(l + z), (1) 
i? 2 = ^P, (2) 

for a flat universe. Here 1 + z is the redshift defined by a/ao = 1/(1 + z). Various linear models [12], [13] have been 
suggested for the EOS in the literatures. We will try to point out in this paper the underlying physics behind these 
linear models and go further to study a more practical leading-orders expansion. 

Taylor expansion is known to be one of the best ways to extract leading contributions from a generic theory with 
the help of a power series expansion of some suitable field variables. The Taylor series is normally convergent due to 
the structure of the expansion coefficients. The power series will, however, converge quickly if the range of variable is 
properly chosen. Hence it is important to choose an appropriate expansion variable for the purpose of our study. 

For example, one may expand the EOS, assuming to be a smooth function for all z, as a power series of the variable 
z around the point z = 0. This will lead to the power series expansion: 

n 

Here the summation of n runs from to 00 and f^ n \z) — d ■ We will not specify clearly the range of summation 
unless it is different from the one we adopted. This convention will be applied to all the summation ranges through 
out this paper for convenience. 

Alternatively, one can also expand the EOS as a power series of the variable w = 1/(1 + z) around the point 
w(z = 0) = 1. The result is 

w(w{z)) = ~i < w - !) = E ~< -H"(t^) (4) 

n n 

Note that if we expand the EOS in power of the variable y = z/(l + z) around y(z = 0) = 0, we will end up with a 
similar power series: 

n 

The Taylor series or the Maclaurin series are normally convergent due to the structure of the expansion coefficients. 
Nonetheless, one would prefer to choose a more appropriate expansion parameter in order to make the series converge 
more rapidly. As a result, leading order terms will be enough to extract the most important physics from the 
underlying theory. Therefore the advantage of the y expansion is that the power series converges rapidly for all range 
of the parameter —1 < y = z/(l + z) < 1, or equivalcntly, — l/2<z<ooas compared to the range \z\ < 1 for the 
power series expansion of z shown in Eq. (3). 

One notes, however, that one should also expand all physical quantities and the field equations to the same order 
of precision we adopted for the EOS expansion. Higher order contributions will not be reliable unless one can show 
that the higher order terms does not affect the physics very much. For example, the liner order is good enough for the 
expansion of the EOS modelling SUGRA model [13]. This is because that the linear term fits the predicted EOS for 
SUGRA model to a very high precision. One readily realizes, however, that it is not easy to track the series expansion 
order by order due to the form of the Eq. (1) for the EOS for the y expansion. This is because that after performing 
the integration, one needs to pay attention to the distorted integration result. Indeed, one needs to write y = 1 — w in 
order to perform the integration involving d\nw. Indeed, one will need to recombine the result back to a power series 
of y. The trouble is that the lower order terms could sometimes hide as the higher order terms in w. It may not be 
easy to track clearly the lower order y expansion when we perform the expansion, for example, due to the exponential 
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factor in the Eq. (1). Therefore, one finds that the most natural way to expand the EOS is to expand it in terms of 
the variable x = ln(l + z) around the point x(z = 0) = 0. Indeed, this power series can be shown to be: 

n 

And this power series converge very rapidly in the range —1 < ln(l + z) < 1, or cquivalently, —0.63 ~ 1/e — 1 < z < 
e — 1 ~ 1.72. Here e ~ 2.72 is the natural factor. Note that this limit happens to agree with the proposed scope of the 
SNAP mission. We will study these different power series expansion for the EOS and its applications in the following 
sections. 



III. POWER SERIES OF Z/(l + Z) 



We will show in details how to extract the leading terms in the y-expansion with y = (j^) around the point 



y(z = 0) = 0. One can write the expansion coefficient as uo n 
expansion of the EOS becomes 



(!+^)<">fa=0) 



l + w(l,(s)) = 5> n y n . 



such that the power series for the 

(7) 



Hence the Eq. (1) can be shown to be 



d\np(y) =3(l+w) 



dy 



= d 



J2 J2 



y 



n+k+l 



n k 



n + k+l 



(8) 



Note that we are now expanding with respect to the smooth function (1 + oj), instead of u), as a Maclaurin series for 
convenience. Therefore, one can integrate above equation to obtain 



P(y) = Pa exp[ 3 Y] - J J -T—, 1 = Po X (y) = Pa Y] X nV n - 



(9) 



n k 



Note that one needs to expand the function p as a power series of y too in order to extract the approximated solution 
with appropriate order. The expansion coefficient X n is defined as X n = X^ n \y = 0)/n\. Here, the superscript in 
X' denotes differentiation with respect to the argument y of the function X(y). One can show that X' = XY with 
Y = 3 J2 n J2 k uj n y n+k . In addition, one can show that 



n k 



(10) 



Hence one has 



Y^(y = 0) = m J2"n- 



n=0 



One can also show, for example, that 



X" = X(Y 2 + Y') 

X'" = X(Y 3 + 3YY' + Y") 

A (4) = X{Y A + 6Y 2 Y' + 3(y') 2 + WY" + Y'"). 



(11) 



(12) 
(13) 
(14) 



This series does not appear to have a more compact close form for the multiple differentiation with respect to y. One 
can, however, put the equations as a more compact format: 



= X[Y + — ] l Y. 
dy 



(15) 
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It appears, however, that one needs to do it manually even it is straightforward. We will only list the leading terms 
as this is already suitable expansion for our purpose at this moment. 
Hence one has 

Xo = 1, (16) 
X! = 3uj , (17) 

X 2 = i[9wg + 3o;o + 3wi], (18) 

X 3 = i[27wg + 27wo(w +wi) + 6(w + wi +w 2 )], (19) 

X 4 = ^[81^o + 162^(w + wi) + 27(wo + t^) 2 + 72w (w + wi + w 2 ) + 18(w + wi + w 2 + w 3 )j. (20) 
Therefore, one can expand the final expression for the energy density p accordingly. Indeed, the result is 
p = p {l + 3co y + ^[9w 2 + 3lv + 3wi]y 2 + ^[27^ + 27w (^o + wi) + 6(w + wi + u) 2 )]y 3 

+ ^[81^o + 162w 2 (w + wi) + 27(wo + ^i) 2 + 72w (^o + wi + w 2 ) + 18(w + tJi + uj 2 + w 3 )]y 4 } (21) 

to the order of y 4 . Note that one keep the order of precision to y A in computing the energy density p even we are 
expanding the EOS only to the order of y 3 . This is due to the special structure in the energy momentum conservation 
law 1. In addition, one can show that the Hubble parameter H = HqX 1 ! 2 with Ho = \J%TtGpo/2>. And the expansion 
for X 1 ! 2 can be obtained by replacing all uj n with uj n /2 in writing the expansion for X . Therefore one has 

3 1 1 

H = H {1 + -uj a y + -[9w 2 + 6(lu + ^i)]y 2 + 7zi 27u o + 54w (w + wi) + 24(w + wi + cj 2 )]y 3 
Z o 4o 

+ y^[162wq + 162w 2 (w + wi) + 54(w + ^i) 2 + 144w (w + u>i + w 2 ) + 72(w + + w 2 + w 3 )]y 4 }. (22) 
Note also that one can also compute the conformal time according to the expression: 

x-i 



Knowing that 1/(1 — y) 2 = J2 n ( n + 1)2/™' onc can show that 

J n 

Therefore, one can easily compute the expansion of 77 in a straightforward manner. 



(23) 



(24) 



IV. POWER SERIES OF ln(l + Z) 



We will show in details how to extract the leading terms in the ^-expansion with x = ln(l + z) around the point 

(l+q;) ( ">(a;=0) 



x(z = 0) = 0. One can write the expansion coefficient as u n 
expansion of the EOS becomes 



l + u(x(z)) = ^oj n x n . 



such that the power series for the 

(25) 



Note that we use the same notation for u n in different parameterizations for convenience. Hence the Eq. (1) can be 
shown to be 



dlnp(x) = 3(1 + ui)dx = d 



n+l 



n+1 



(26) 
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Note that we are now expanding the physical quantities with respect to the function (1+w) instead of u> for convenience. 
Therefore, one can integrate above equation to obtain 



n+l 

p(x) = po exp[ 3 ^ -^-^ } ee Po X(x) = p Q J] X n x n . (27) 

n n 

Note that one needs to expand the function p as a power series of x too in order to extract the approximated solution 
with appropriate order. The expansion coefficient X n is defined as X n = X^ n \x = 0)/n!. One can show that 
X' = XY with Y = 3 J2 n = 3 ( : + w )- Therefore, one has 

YW(y = 0) = 3(l\)uj l . (28) 

One can also show, for example, that 

X" = X{Y 2 + Y') (29) 
X'" = X(Y 3 + 3YY' + Y") (30) 

x^ = x(y 4 + 6y 2 y' + 3(y') 2 + Try" + Y'"). (31) 

In addition, one can show that 

n fe v ' 

Therefore, one has 

X®= ^ 3X„w fc (i-l)!. (33) 

n+fc=/-l>0 

Hence one obtains the recurrence relation for the expansion coefficients of X n : 

3 l-i 

x, = y5353x„ Wfc . (34) 

As a result, one has, for example, 

X = 1, (35) 
Xi = 3w , (36) 

A- 2 = ^[9wg + 3wi], (37) 



A 3 = - [27wg + 27uwi + 6lo 2 ] , (38) 



X 4 = t^; [81^o + 162wqWi + 27w 2 + 72w w 2 + 18w 3 ]. (39) 
Therefore, one can expand the final expression for the energy density p accordingly. Indeed, one has 

p = p {l + 3uj y + ^[9lu 2 + Suj^y 2 + ^[27w% + 27u ^i + 6w 2 ]y 3 
+ ^i 81uj o + 162wgwi + 27ujj + 72uj uj2 + 18w 3 ]y 4 }. (40) 



In addition, one can show that the Hubble parameter H = HqX 1 / 2 with H — y / 87rGpo/3. And the expansion for 
A 1 / 2 can be obtained by replacing all u n with w n /2 in writing the expansion for A. The result is 

3 1 1 

H = H {1 + -cj y + -[9uj 2 + 6wi]y 2 + —[27^ + 54w ^i + 24w 2 ]y 3 

+ -!- [81wJ + 324wjjwi + 108w 2 + 288^2 + 144w 3 ]y 4 }. (41) 

oo4 
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Note also that one can also compute the conformal time according to the expression: 

% = («-U/^. (42) 
Knowing that 1/(1 — y) 2 = J2 n ( n + 1)2/™) one can show that 

H oV = I dyX-^Y,{n+l)y n . (43) 

J n 

Therefore, one can easily compute the expansion of rj in a straightforward manner. 

V. CONCLUSION 

The proposed Supernova/ Acceleration Probe (SNAP) will carry out observations aiming to determine the compo- 
nents and equations of state of the energy density, providing insights into the cosmological model, the nature of the 
accelerating dark energy, and potentially clues to fundamental high energy physics theories and gravitation. As a 
result, we are motivated by a nice result shown in Rcf. [13] to study the physics underlying the model presented there. 

A detailed discussion on the choices of the expansion parameters of the Maclaurin series for the equation of states 
of a perfect fluid is presented in this paper accordingly. For example, the Maclaurin series of the EOS is expanded 
as power series of the variables y = z/(l + z) and x = ln(l + z) respectively. We also show how to obtain the power 
series for the energy density function, the Hubble parameter and related physical quantities of interest. The method 
presented here will have significant application in the precision distance-rcdshift observations end to map out the 
recent expansion history of the universe, including the present acceleration and the transition to matter dominated 
deceleration. We also show that the nice recent result in Ref. [13] is in fact a linear approximation to the full Maclaurin 
series as a power series of the parameter y = z/(l + z). 

Since we can power expand all smooth EOS into a convergent power series, it is more practical for the future probe 
to measure local derivatives of the EOS. Note that u>o represents the value of 1 + u at time to or z = 0. In addition, 
lo\ is the slope of the curve u) vs y or u> vs x at z = 0. Therefore, knowing the values of the expansion coefficients 
will enable us to plot the entire EOS through out the range of convergence. One may need to use different expansion 
series depending on the convergent speed of the power series. For example, it appears that the leading order term in 
the y expansion is good enough to obtain a nice result for the SUGRA prediction. This is because the leading term 
is close enough to the theoretical prediction [13]. Nonetheless, one expects that a few leading terms in the Maclaurin 
series will be able to offer close enough result for us in the future. 

In addition, the result shown here is independent of the choice of the time to- The local measurement of the 
expansion coefficients can be extended to compare the expansion coefficient at any time. One probably should trust 
less on the results other than the local predictions. 
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